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ABSTRACT 

The benchmark problem 3 in Category 3 of the third 
Computational Aero- Acoustics (CAA) Workshop sponsored by 
NASA Glenn Research Center is solved using the space-time con- 
servation element and solution element (CE/SE) method. This 
problem concerns the unsteady response of a rectilinear swept 
cascade to an incident gust. The acoustic field generated bv the 
interaction of the gust with swept flat plates in the cascade is 
computed by solving the 3D nonlinear Euler equations using the 
space-time CE/SE method. A parallel version of the 3D CE/SE 
Euler solver is employed to obtain numerical solutions for several 
sweep angles. Numerical solutions are presented and compared 
with the analytical solutions. 

INTRODUCTION 

Fan noise is one of the dominant sources of engine 
noise on both takeoff and landing approach for turbofan 
engines with high bypass ratio. Thus the development of 
fan noise reduction technology is very important. Noise 
source control is one of the strategies to reduce the fan 
noise. One of the significant fan noise sources is the rotor- 
stator interaction noise. The fan (rotor) blades generate vis- 
cous wakes that appear as the periodic velocity fluctuations 
that are convected by the mean flow to interact with the 
downstream outlet guide vanes (stator), producing acous- 
tic waves that propagate upstream and downstream. The 
benchmark problem 3 in Category 3 of the third CAA Work- 
shop is a model problem of rotor-stator interaction noise. 
The gust represents the rotor viscous wake and the cascade 
of swept flat plates represents the stator. The effect of the 
sweep angle of the flat plate on the upstream radiated noise 
is investigated. 


This benchmark problem is solved numerically by solv- 
ing the unsteady 3D Euler equations using the space- 
time conservation element and solution element (CE/SE) 
method. The CE/SE method is an innovative numerical 
method for solving conservation laws. This method is dis- 
tinguished from other methods by its requirement that the 
fluxes be conserved in both space and time. Simplicity, gen- 
erality and accuracy were the goals in the development of 
this method. Its salient properties are summarized briefly as 
follows. First, both local and global flux conservations are 
enforced in space and time instead of in space only. Second, 
all the dependent variables and their spatial derivatives are 
considered as individual unknowns to be solved for simulta- 
neously at each grid point. Third, every CE/SE scheme is 
based upon a non-dissipative scheme with addition of fully 
controllable numerical dissipation. This results in very low 
numerical dissipation. Fourth, the flux-based specification 
of the CE/SE schemes gives rise in a natural fashion to a 
simple yet generally effective non-reflecting boundary con- 
dition which is an important issue in CAA. Some practical 
advantages of the CE/SE method over the high order fi- 
nite difference methods are summarized as follows: 1) it 
can handle both shocks and acoustic waves concurrently, 2) 
it can use both structured and unstructured meshes in one 
single algorithm to handle complex geometries, 3) it can 
avoid singular points without using any special treatment, 
4) it has the most compact stencil, this leads to efficient 
parallel computing and easy implementation of boundary 
conditions, and 5) it can capture shocks accurately without 
using Riemann solvers or dimensional-splitting techniques. 
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A detailed description of this method and accompany- 
ing analysis can be found in [1-5]. Applications of this 
method to CAA problems reveal that, on mesh sizes used 
in practice, the result is comparable to that obtained using 
high order compact difference schemes even though the cur- 
rent solver is only 2nd-order accurate [6-11]. Especially for 
the CAA benchmark problems involving both shocks and 
acoustic waves, the CE/SE method can produce accurate 
solutions with much less effort than the high order finite 
difference methods [10,11]. 

In t his paper, a parallel version of the 3D CE/SE Euler 
solver is used to compute the acoustic field generated by 
the interaction of a gust with a rectilinear swept cascade, 
which is described as the benchmark problem in Category 3 
of the third CAA Workshop [12]. Different sweep angles of 
flat plates in the range of [0°,30°] are studied to determine 
the upstream-radiated noise as a function of sweep angle 
at the specified frequency. The variation of RMS pressure 
at an upstream point with the sweep angle is plotted and 
compared with the analytical solution. 

In the following, the description of the benchmark prob- 
lem is given first, which is followed by the boundary and ini- 
tial conditions used in the calculation, and numerical results 
and discussion. 


GUST-CASCADE PROBLEM 

Consider the rectilinear cascade of swept flat plates 
shown in Fig. 1. The x-axis is aligned with the chord of 
the airfoil, the y - axis is perpendicular to it and the £-axis 
is normal to the bounding walls. The normal distance be- 
tween w r alls is 2.6c, where c is the chord length of the flat 
plate. The gap-to-chord ratio, sfc, is 1.0. The sweep angle 
of the flat plate, denoted by 0. varies from 0° to 30°. 

The mean flow is assumed to be uniform and aligned 
with the .r-axis. The mean flow variables are inflow ve- 
locity, Co, static density, po, and static pressure, po- The 
inflow Mach number, A/o, is 0.5. Flow variables are non- 
dimensionalized by using ao(the speed of sound) as the ve- 
locity scale, c as the length scale, c/ao as the time scale, po 
as the density scale and po< 2 o as the pressure scale. Thus 
the mean flow is described in dimensionless variables as 

p — 1 , u = 0.5, v = 0, w — 0, p=l. 0/1.4 (1) 

The incident gust carried bv the mean flow has ar, y , 
and z velocity components given by 

u — — {vak y /k x ) cos (k x x -rk y y~rk z z — cJt) (2) 


and 

v* = va cos(k x x + k y y — k z z - art), w f = 0 (3) 

respectively, where vq — 0.0001 and k x = 5.5 , — i r,fc c = 
0, and lj = Mok x , respectively. The corresponding period 
of the gust wave is T — 2i t/oj. It is assumed that the gust 
is '‘frozen” and convected by the uniform mean flow. Thus 
the gust satisfies the linearized Euler equations. This im- 
plies that p' = p' = 0. The gust is continually convected 
from the inlet by the mean flow and will interact with the 
swept flat plates to generate acoustic waves. The whole flow 
field, including the acoustics, gust wave and mean flow, is 
simulated by solving the nonlinear Euler equations. 

BOUNDARY AND INITIAL CONDITIONS 

Having the gust superposed on the mean flow, the flow 
properties at the inlet are described as: 


(«i )j=P (4) 

M” = p{u -r u') (5) 

(«••?)" = P{v + <’') (6) 

(“ 4 )” = p(w + w') (7) 


(« 8 )? = 7TT - f I [(“ + u'f + (v + v'f + {w + «/) 2 ](8) 

The inlet values of ( u m x )^, (umy )™ 5 and (u mz )j are ob- 
tained by evaluating derivatives of u m at (j, n). This is one 
method for imposing the non-reflecting BC in the CE/SE 
method [6]. An alternative non-reflecting scheme is used at 
the outlet, as follows 

(u m y; +l = Kntfty 2 , (u mx y; +l = 0 

(n my )] +1 = (Umy)”±t /2 , (Um:)] +1 = (9) 

This way of imposing the non- reflecting boundary condi- 
tions allows the flux to “stream” out of the spatial domain 
smoothly with minimal reflection. The periodic boundary 
condition is imposed on planes y = —0.5 and y = 1.5 and 
the reflecting boundary conditions are used on the flat plate 
surfaces and top and bottom walls in the ^-coordinate direc- 
tion. No grid points are located at the flat plate leading and 
trailing edges, in order to avoid the singular flow behavior. 

At t = 0, the time-marching variables (u m ) 7 j in the 
whole domain are defined using Eqs. (4)-(8) with t — 0, 
and (Umar)” , and (u mz )j are obtained by evaluating 

derivatives of u m at (j, n). 
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NUMERICAL RESULTS AND DISCUSSION 

A structured 291x21x27 hexahedral grid is generated 
by using an algebraic transformation in the computational 
domain of —6 < x < 8.5, —0.5 < y < 1.5. and 0 < z < 2.G 
first. A slice of the grid in x-z planes is shown in Fig. 2. 
Each hexahedron is cut into six tetrahedrons in the way 
shown in Fig. 3, which results in 904800 (290x20x26x0) 
tetrahedral cells in the computation. The parallel version of 
the 3D Euler solver is used with a — 0 and e — 0.5, where a 
and e are the parameters connected with the specification 
of the numerical dissipation in the solver [4]. A detailed 
description of the 3D Euler solver is referred to in [4] and the 
parallelization part is given in [5]. For 9 < 25°, A t = T j 168, 
while At = T/ 200 for 9 = 30°. 

The sweep angle of flat plates is varied to determine its 
effect on the upstream-radiated noise level. Sweep angles 
in the range of [0°,30°] are studied here. In all figures, the 
acoustic pressure p f non-dimensionalized bv vq is plotted. 
The time history of the acoustic pressure at an upstream 
point located at around (-5, 0,1.3) is shown in Fig. 4 for 
sweep angles in the range of [2°, 30°]. It should be pointed 
out that numerical instability arises at leading edge with 
the current meshing strategy shown in Fig. 3 for sweep an- 
gles that are less than 5°. No converged solution is obtained 
for 6 = 0°. Some instabilities are observed for 9 = 2°. This 
instability can be overcome by using a different meshing 
strategy, such as cutting each hexahedron into 24 tetrahe- 
drons through diagonal lines, which is not presented here. It 
can be seen that the solution is fully converged for 9 > 5° by 
t = 46T = 105.1. And the absolute amplitude of the acous- 
tic pressure at the upstream decreases when the sw'eep angle 
increases. 

The numerical result of the RMS pressure at the up- 
stream point needs to be expressed in dB using the value 
for 9 — 0° as the reference. However, the computed result of 
the RMS pressure for 6 = 0° is not available. Thus the value 
of the RMS pressure for <9 = 5° is used as the reference with 
a correction such that the value of the computed dB pres- 
sure is identical to the analytical solution for 6 = 5°. The 
so-obtained data is plotted in Fig. 5 along with the ana- 
lytical solution [13] that shows the upstream-radiated noise 
level is reduced substantially by increasing the sweep angle, 
and the noise field is cut-off at the critical sweep angle. A 
reasonable agreement between numerical and analytical so- 
lutions is observed for smaller sweep angles (< 15°), while 
there are significant discrepancies for larger sweep angles 
(> 20°). One possible explanation is that the effect of the 
3D corners formed by the upper (lower) wall and the flat 
plates becomes significant for larger sweep angles. However, 


the 3D corner effect is not fully taken into account in the 
analytical solution [14]. The numerical error caused by the 
reflection at the outlet is investigated by running the code 
in a smaller domain. A similar grid of 201x21x27 in a do- 
main of -6 < x < 4, -0.5 < x < 1.5, and 0 < z < 2.6 
is used with the identical boundary conditions. The cor- 
responding data is plotted for 9 - 5°, 10°, and 15° in Fig. 
6. As expected, a slightly better result is produced by us- 
ing the larger domain. The solutions obtained using two 
different domains are very close, which means the bound- 
ary condition has small effect on the accuracy of numerical 
solutions. 

Further, the acoustic pressure contours on the upper 
and lower surfaces of the flat plate are plotted in Fig. 7 
for several sweep angles to show the different wave patterns 
near the noise source. The pressure contours on x-y and 
y-z planes and pressure distribution at some constant lines 
are shown in Figs. 8-9 for two sweep angles. It can be seen 
that the variation of the acoustic field in the z- direction is 
amplified by increasing the sweep angle and the solutions 
show the periodicity at the periodic boundary. 

The parallel 3D CE/SE code was run on Origin2000 
clusters. For a mesh of 904800 cells, it takes around 5 hrs 
of wall-clock time to reach t = 46T (15456 time iterations) 
using 31 CPUs. 


CONCLUSION 

The problem 3 in Category 3 of the third CAA Work- 
shop has been solved by using the 3D parallel CE/SE Euler 
solver. Numerical results are presented for several sweep 
angles. The result of the RMS pressure at an upstream 
point is compared wdth the analytical solution, showing a 
reasonable agreement for smaller sweep angles (< 15°). 
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Figure 1. GEOMETRY OF SWEPT RECTILINEAR CASCADE IN A FINITE 
DOMAIN. 



Figure 2. A SLICE OF THE HEXAHEDRAL MESH ON THE X-Z PLANE. 


E 



Figure 3. SIX TETRAHEDRONS IN EACH HEXAHEDRON. 
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(c) slices at i- constant lines 





Figure 8. ACOUSTIC PRESSURE CONTOURS AND DISTRIBUTION ON X - Y AND Y - Z PLANES FOR 0 = 5 



Figure 9. ACOUSTIC PRESSURE CONTOURS AND DISTRIBUTION ON X - Y AND Y - Z PLANES FOR 9 = 25 
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